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Abstract 

Atmospherically-deposited nitrogen (N) can stimulate complex soil N metabolisms 
and accumulations over time. Whether long-term (decadal) N deposition effects on 
soil N transformations and functional microbes differ from the short-term (annual) 
effects have rarely been assessed. Here we conducted a laboratory PN tracing study 
with soil samples from a short-term (one year) N addition site and a long-term (12 
years) site in a subtropical forest. The effects of simulated N deposition on soil N20 
emissions, N transformation rates and microbial nitrifying and denitrifying genes 
were determined. Our results showed that: (1) long-term N addition did not change 
soil №0 fluxes significantly in comparison to the short-term М addition. 
Denitrification, heterotrophic nitrification and autotrophic nitrification contributed 
53%, 23% and 18% to total №0 emissions, respectively. (2) Autotrophic nitrification 
was the dominant N transformation process, except for the high-N treatment at the 
long-term site. The magnitude of soil N transformation rates was significantly 
different among N addition treatments but not between short- and long-term N 
addition sites. However, long-term N addition changed the responses of specific N 
transformation rates to N addition markedly, especially for the rates of nitrification, 
organic N mineralization to МНД”, NO3 immobilization and dissimilatory NOx 
reduction to МН4 (ОМКА). (3) Responses of ammonia oxidizing archaea and bacteria 
(AOA and АОВ) were stronger than denitrifying N2O-producers (nirK) and 
denitrifying N2O-reducers (nosZ) at the long-term site compared to the short-term site. 
(4) The close correlations among №0 flux, functional genes and soil properties 
observed at the short-term site was weakened at the long-term site, posing a decreased 
risk for N losses in the acid subtropical forest soils. There is evidence for an 


adaptation of functional microbial communities to the prevailing soil conditions and 
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in response to long-term natural and anthropogenic N depositions. 
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1. Introduction 


Nitrous oxide (№0) is an important greenhouse gas with stronger radiative 
forcing than СО» (Matson et al., 2002; Wang et al., 2014a). Terrestrial ecosystems 
contribute approximately 57% of the global N2O emissions (Werner et al., 2007). In 
particular, tropical and subtropical forests act as the largest natural terrestrial source of 
N20 globally (Werner et al., 2007; Cheng et al., 2014). In recent years, tropical and 
subtropical regions have been experiencing severe atmospheric N deposition due to 
the enhanced atmospheric reactive М (N;) (Kanakidou et al., 2016) mainly come from 
anthropogenic N fertilization (MacDonald et al., 2002; Herrmann et al., 2005; Du et 
al., 2014; Zhu et al., 2015), biomass burning and fossil fuel combustion (Boy et al., 
2008; Chen et al., 2010). The deposited inorganic and organic N are subsequently 
involved in soil internal М turnovers (Avrahami et al., 2002; Bai et al., 2014; Deppe et 
al., 2017). Forest soil N20 emissions are therefore likely affected by the increased 
atmospheric N deposition (Galloway et al., 2004; Galloway et al., 2008; Vet et al., 
2014). However, the relationship between increased N deposition and soil N 
transformation/N2O emission in subtropical and tropical forests is poorly understood 
due to the complex interactions among various influencing factors such as climate, 
forest type, soil property, microbial community and land-use history (Chatskikh et al., 
2005; Alm et al., 2007; Werner et al., 2007; Zhang et al., 2008; Chen et al., 2014). 

Although abiotic processes such as chemodenitrification and chemical 
decomposition of hydroxylamine may contribute to soil N20 production (Williams et 
al., 1992; Bremner, 1997; Butterbach-Bahl et al., 2013), microbial nitrification and 
denitrification contribute over 70% of the global N2O emissions (Avrahami et al., 


2002; Braker and Conrad, 2011; Butterbach-Bahl et al., 2013). Different groups of 


soil microorganisms are involved in Ше processes of soil М transformation that 
produce N20 (Levy-Booth et al., 2014). For example, ammonia (NH3) oxidization, 
the first step of autotrophic nitrification (Омн4) that catalyzed by the amoA-encoded 
ammonia monooxygenase of ammonia oxidizing archaea and bacteria (AOA and 
AOB), is an important rate-limiting step for N2O production (Szukics et al., 2010; 
Isobe et al., 2012; Long et al., 2012; Wertz et al., 2012; Wang et al., 2014b; Tang et 
al., 2016). Denitrification, i.e. reduction of nitrate (NO3°) to №0 and further N20 to 
№ that catalyzed by several enzymes such as narG-encoded nitrate reductase, 
nirK/nirS-encoded nitrite reductase, norB-encoded nitric oxide reductase and nosZ- 
encoded nitrous oxide reductase, is also essential for the formation and consumption 
of N20 in soils (Freedman et al., 2013; Jones et al., 2013; Levy-Booth et al., 2014; 
Orellana et al., 2014; Yu et al., 2014; Domeignoz-Horta et al., 2017). Quantification 
of soil nitrifiers and denitrifiers, especially nitrifying N2O-producers (АОА-атоА, 
AOB-amoA), denitrifying N2O-producers (nirK/nirS) and denitrifying N2O-reducers 
(nosZ and the newly identified N2O-reducing clade), may facilitate the understanding 
of N20 productions and consumptions in soils (Isobe et al., 2012; Yu et al., 2014; 
Zhang et al., 2014; Domeignoz-Horta et al., 2015; Tang et al., 2016; Domeignoz- 
Horta et al., 2017; Hallin et al., 2017). Moreover, a range of environmental factors, 
such as soil moisture, redox potential, pH and nutrient availability, have been proved 
to be able to affect soil nitrifiers and denitrifiers (Bárta et al., 2010; Hu et al., 2013; 
Zhang et al., 2014; Zhang et al., 2015; Tang et al., 2016). Although a number of 
studies have been devoted to understanding the relationships among soil properties, 
functional microbes and №О fluxes (Fang et al., 2008a; Yu et al., 2014; Zhang et al., 
2014; Faeflen et al., 2016; Chen et al., 2017), our knowledge for the responses of 


microbial-regulated №0 emission to elevated М deposition in subtropical forest soils 


is still limited increasing backgrounds is still limited (Avrahami et al., 2002; Zhong et 
al., 2015; Cui et al., 2016; Gao et al., 2016). 

Soil №0 production and emission occur in diverse soil М pools via different М 
transformation processes, including Ше №0 originates from the nitrate nitrogen (МОз 
-N) pool via denitrification, from the ammonium nitrogen (ХН4 -Х) pool ма 
autotrophic nitrification and from the organic nitrogen (Norg) pool via heterotrophic 
nitrification (Müller et al., 1998; Zhang et al., 2011). Quantifying the contributions of 
these №0 production pathways is essential to understanding the fate of 
atmospherically-deposited inorganic N and organic N (Fang et al., 2008b). Since the 
production and consumption of NH4* and NO: are also tightly associated with the 
mineralization of organic N to NH4* and МОз, the immobilization of NH4* and МОз 
to organic М, and the dissimilatory NO3 reduction to МНА (ОМКА), quantifying the 
contributions of these processes is also crucial to uncovering the fate of 
atmospherically deposited N. In the last decade, progress has been made via advanced 
microbial and stable isotope methods to more accurately understand internal N 
transformations and underlying soil microbial drivers. The development of РМ tracing 
techniques allows the quantification of specific N processes by differentially labelling 
soil N pools (e.g. ^NH4* ог Оз) (Hart et al., 1994). Moreover, the 
atmospherically-deposited N will stimulate complex soil N metabolisms and 
accumulations over time (Fang et al., 2009; Wang et al., 2015; Deppe et al., 2017). 
The effect of N deposition over long periods of time will likely induce rate and 
microbial responses that are different from short-term effects. Some previous studies 
have addressed the temporal effects of N deposition on forest soil N transformations 
but mainly focused on seasonal comparisons (Breuer et al., 2000; Kiese et al., 2003; 


Zhu et al., 2013a); less attention has been paid to the differentiation of short-term (a 


couple of years) vs. long-term (over a decade) N deposition effects (Isobe et al., 2012; 
Gao et al., 2016; Gurmesa et al., 2016). 

Subtropical forests in south China are suffering from high atmospheric N 
deposition (Mo et al., 2006; Zhu et al., 2013a; Zhu et al., 2015) with hitherto 
unknown effects on ecological processes. The Dinghushan Long-term Nitrogen 
Research (DHSLTNR) study (Lu et al., 2008; Liu et al., 2011) is one of the longest- 
running N deposition projects in China and has received 12 years of simulated N 
addition. This provides a unique research facility to conduct time-scale comparisons. 
To compare the long- vs. short-term effects of N deposition on soil N transformations 
and associated functional microbes, we further established a short-term N addition site 
that is adjacent to the long-term site. The similar hydrothermal conditions at the two 
sites allow the comparison of N2O emission, N transformation as well as gene 
regulation differences with time at different N deposition levels. Our specific 
objectives of the study were to: 1) quantify the relative contribution of N20 
production pathways in the subtropical forest soils, 2) characterize the responsive 
behaviors of various N transformation processes to simulated N deposition, 3) assess 
the role of functional microbes in regulating N20 emission and the fate of deposited 


N, and 4) compared the long- vs. short-term N addition effects. 


2. Materials and methods 


2.1 Site description and sample collection 


The simulated short-term and long-term N deposition sites are located in the 


Dinghushan Biosphere Reserve (DHSBR) (23?10' М, 112?10' E) which belongs to 
Guangdong province of south China. The DHSBR covers 1155 ha, and the climate 
there is characterized as humid monsoon with 1927 mm mean annual precipitation. 
About 75% of the precipitation occurs from March to August with a distinct wet-dry 
seasonality. The mean annual temperature of this zone is 21.0 °C, with a minimum 
and maximum monthly average value of 12.6 °C in January and 28.0 °C in July (Mo 
et al., 2006). In addition, the reserve has experienced naturally high rates of 
atmospheric N deposition (20-50 kg М Па а!) with total wet N deposition reaching 
34.4 kg N Ва а! compared to the average value of 29 kg N ha’! a`! (as inorganic N in 
bulk precipitation) in the 1990s (Fang et al., 2008b; Lu et al., 2013; Gurmesa et al., 
2016). The two sites studied in this article are located in a monsoon evergreen 
broadleaf forest which has been well protected from human activities in the past 400 
years (Mo et al., 2006). The soil is Ferrasols based on the World Reference Base for 
Soil Resources (WRB) classification system (FAO, 1998; Chen et al., 2005) with the 
depth ranging from 30 to 70 cm (Mo et al., 2006; Fang et al., 2009). Major tree 
species in this monsoon forest include Castanopsis chinensis, Schima superba, 
Gironniera subaequalis, Syzygium acuminatissimum, and Aporusa yunnanensis (Mo 
et al., 2006). 

The field N addition experiment at the long-term site (L site) was initiated from 
July 2003 by spraying ammonium nitrate (NH4NO») solution monthly (Mo et al., 
2006). Four N-addition gradient treatments were set up with applications of 0, 50, 100 
and 150 kg N Ба a`! as NH4NOs. Hereafter, the four treatments are abbreviated as: 
control (LC), low N (LL), medium N (LM) and high N (LH), respectively. Three 
replicate plots (each 10 m x 20 m) were established for each N treatment (in total 12 


plots). A buffer strip of about 10 m was left between adjacent plots. For each N- 


addition plot, Ше required NH4NO3 amount to Ше target application was dissolved in 
20 L water and sprayed evenly below the canopy using a backpack sprayer (1.е. 
application rate of 0.1 L m? was low to avoid liquid application effects) (Mo et al., 
2006; Lu et al., 2013). The short-term site (S site) received N depositions from 
September 2014 onwards at 0, 35, 70 and 105 kg М ha’! a! for the control (SC), low N 
(SL), medium N (SM) and high N (SH), respectively. Three replicates (each 15 m x 
15 m, with a 10 m wide buffer strip between every adjacent two plots) under each N 
treatment formed 12 plots in total. МНАМОз fertilizer was dissolved in 30 L water and 
sprayed evenly to N-plots (i.e. application rate of 0.1 L m?). At the time of soil 
sampling in July 2015, the long-term site and the short-term site had received 
simulated N deposition for approximately 12 years and one year, respectively. 
Twenty-four soil samples (2 sites х 4 М treatments x 3 replicates) were collected 
from the surface soil (0 - 20 cm) after litter removal in July 2015. Six soil cores were 
randomly selected in each replicate plot and mixed into a composite sample. АП soil 
samples were sieved (2 mm mesh) and split into three subsamples: one subsample 
was stored at 4 ?C for the incubation studies, the second subsample was stored at - 
80 °С for functional gene analyses, and the third subsample was directly used or air- 
dried for soil physiochemical property determination, including soil water content 
(WC), water holding capacity (WHC), soil pH, ammonium nitrogen (NH4'-N), nitrate 
nitrogen (NO;-N), total organic carbon (TOC), total nitrogen (total-N), microbial 


biomass carbon (MBC) and nitrogen (MBN). 


2.2 ÜN-tracing experiment 


A paired PN tracing experiment was conducted for the collected soil samples. 


There were two PN treatments: РМНАМОЗ and NH4!°NOs, both at 10 atom% ^N 
enrichment. For all soil samples collected from the two sites, 240, 250 ml flacks (2 
'SN treatments x 2 sites x 4 N treatments x 3 replicates x 5 sampling points) were 
prepared. Each flask contained fresh soil (equivalent to 20 g dry weight) and was 
covered by a sealing membrane with three pin-holes to allow gas exchange but avoid 
moisture losses. Soils were incubated in the dark at 20 °C and 40% WHC for an 
overnight pre-incubation. After that, 2 ml РМНАМОз or МНАРМОз solution was 
evenly added to the soil at an amount of 20 mg N kg"! soil (oven-dried weight), and 
the labelled soils were incubated in the dark at 20 °C and 60% WHC for 336 h (14 d). 
Gas samples were taken from the headspace of the flasks at 0.5 h (0 d), 24 В (1 d), 72 
h (3 d), 168 h (7 d) and 336 h (14 d) after label addition. The incubation and gas 
sampling were carried out according to Zhang et al (2011). Finally, two 30-ml gases 
were taken from the headspace of each flask with a gastight syringe (30 mL, Ping An, 
China) and transferred into pre-evacuated vials (12 mL, Jing Yang, China). One gas 
sample was used to determine №0 concentration while the other sample was used for 
the analysis of ^N enrichment in №0. After gas collection, soil samples were 
extracted with 2M КС solution at a ratio of 1: 5 (soil: extractant). After KCl 
extraction soils were washed with distilled water to exclude inorganic nitrogen, then 


filtered and dried for the analysis of ^N enrichment in organic М. 


2.3 Determination of soil properties, N20 fluxes and МО pathways, N transformation 


rates and gene abundances 


2.3.1 Soil properties 


Soil NH4'-N was measured using the Alpha-Naphthol Blue- spectrophotometer 
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method after KCI extraction, and NO3-N was measured using Ше Dual Wavelength 
Spectrophotometric method. Soil MBC and MBN were measured using the 
chloroform-fumigation-extraction method (Vance et al., 1987). Soil pH was measured 
using the mixture of soil/water at a ratio of 1:2.5 (т: у) Буа portable pH detector (F- 
71G, LAQUA, HORIBA, Japan). Soil TOC was measured using the H2SO4-K2Cr207- 
heating method. Total-N content was analyzed after H2SOx digestion using both the 
Alpha-Naphthol Blue- spectrophotometer method and the Mo-Sb Anti- 


spectrophotometer method (Liu, 1996). 


2.3.2 N concentrations and М enrichments in inorganic-N, organic-N and N20 

The РМ enrichment of soil inorganic М (NH4* and МОз`) in the KCI extract was 
measured by an isotopic ratiomass spectrometer (IRMS, Europa Scientific Integra, 
Crewe, UK). УМНА? and PNOs' were first separately transferred into NH; and then 
converted into (IPNH4)2SO4 (Wang et al., 2015), while the concentrations of the two 
inorganic-N forms were measured Бу 0.005 М H5SO. titration using a 2% НзВОз 
solution as a color indicator. Briefly, the KCl extract (containing ^NH4* and PNO;) 
was steam-distilled with magnesium oxide (MgO) for converting NH4 to liberated 
NH; and the same extract in the flask was distilled again with Devarda’s alloy for 
transferring МОз to liberated МН; (Feast and Dennis, 1996). The liberated NH3 
produced in the twice distillations was separately trapped using the boric acid solution 
(i.e. a titration for measuring the concentration of NH4* and МОз). The trapped N in 
the flask was then acidified by the Н5О4 solution for converting to (NH4)2SO4. The 
solution after acidification was dried at 80?C for ^N enrichment determination. After 
KCI extraction and filtration, the soils left were washed four times with distilled water 


and subsequently dried at 60 °С for the determination of organic-N concentration and 
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ISN enrichment. For Ше gas samples, N2O concentration was measured using ап 
Agilent 7890A gas chromatograph (Agilent Technologies, Inc, USA), and the ^N 
enrichment in №0 was measured using a mass spectrometer (Finnigan, MAT 253 


mass spectrometer). 


2.3.3 Calculation of N20 fluxes and N20 pathway contributions 

The №0 flux in each sampling point was calculated according to standard 
methods (Yu et al., 2014; Xie et al., 2015) based on the N20 concentration change 
over time during the incubation. Based on the РМ enrichment (atom% excess) in ХО, 
NH4'-N, NO3-N and also in organic-N pools in the paired РМ tracing treatments, we 
calculated the contribution of different N2O emission fractions. We assumed that N20 
originated from NH4'-N, NO3-N and organic М pools ма autotrophic nitrification, 
denitrification and heterotrophic nitrification, respectively. The contributions of these 
three different pathways were calculated according to the following equations 


(Rütting et al., 2010; Zhang et al., 2011): 

аро 7d:a, аа t h-a, (1) 

а+а+ћ =1 (2) 
Where ayo , 0j, d, and а, indicate Ше РМ atom% excess in N20, МОз, NH4*, and 
organic N іп ће '"NH4*-labelled ог '"NO3-labelled treatment during the incubation; 
d, aand h indicate the N2O emission fraction that originates from the NO3-N pool via 
denitrification, from the NH4'-N pool via autotrophic nitrification and from the 
organic М pool via heterotrophic nitrification. By inserting relevant РМ enrichment 


data into the equations above we calculated the values of d, a and h by the back solver 


method in Excel (Microsoft, Inc.). 
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To understand ^N losses from Ше labelled inorganic N pools in form of №0, we 
calculated the percentage of the РМ enrichment (atom% excess) in N20 to the ^N 
enrichment in the NH4*-N pool, i.e. ^N20/PNH4*96 in the '"NHa*-labelled treatment, 
and also the percentage of the ^N enrichment іп N20 to the PN enrichment in the 
NO3-N pool, i.e. ^N20/PNOs % in the "NO3-labelled treatment for each sampling 


point of the incubation. 


2.3.4 Quantification of N transformation rates 

Gross N transformation rates were quantified by а '°N tracing model (Müller et 
al., 2007) (Fig S1). The model considered ten transformation processes related to the 
turnover of organic М and inorganic М in soil: (1) Ммес, mineralization of recalcitrant 
organic N to NH; (2) Гуна-мес, immobilization of МНА to recalcitrant organic М; (3) 
Ммаь, mineralization of labile organic N to МН4; (4) Гунч-маь, immobilization of 
МН. to labile organic М; (5) Омес, oxidation of recalcitrant organic N to NO3 
(heterotrophic nitrification); (6) Гуоз, immobilization of NO: to recalcitrant organic 
N; (7) Охна, oxidation of NH4' to NO; (autotrophic nitrification); (8) Омоз, 
dissimilatory NOs reduction to МНА (ОМКА); (9) Амна, adsorption of NH4' on cation 
exchange sites; (10) нула, release of absorbed МНА. The tracing model calculated 
these N transformation rates based on the concentration and PN enrichment values 
(average + standard deviations) of inorganic М and organic М from the two distinct 
'SN labelled treatments during the experimental duration. Finally, N transformation 
rates of all the eight N addition treatments at the two sites were calculated over the 14 
d period and expressed in units of mg М kg" soil day’. 

Based on the results from the ^N tracing model we calculated net NH4* 


production rates (МЕТхнарпоа.) and net NO3 production rates (МЕ Two3prod.): 
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NET на prod. — M Nrec T M Nlab T Кунда +D NO3 ` I NH4-Nrec — I NH4-Nlab Avia (3) 


NET, NO3prod. — Oy, + Ong, —1 №3 — Dyos (4) 


2.3.5 Soil DNA extraction and functional gene abundance quantification 

DNA was extracted from 0.25 g fresh soil (before incubation) using PowerSoil° 
DNA Isolation kit (ANBIOSCI TECH LTD) following Ше manufacturer’s 
instructions. The extracted DNA solution was checked for quantity and quality using a 
Nanodrop spectrophotometer (NanoDrop Technologies Inc., Wilmington, DE). 
Meanwhile, AOB-amoA and AOA-amoA were chosen to represent nitrifying 
community groups, nirK and nosZ were chosen to represent denitrifying community 
groups. The primers of four genes are listed in Table S1 based on previous studies 
(Rich et al., 2003; Levy-Booth et al., 2014). 

Quantitative PCR was carried out with an ABI 7500 CFX96 Optical Real-Time 
Detection System (Bio-Rad Laboratories, Inc., Hercules, CA) by using SYBR green 
kits. The 20 uL reaction mixture contains 10 uL SYBR Premix Ex Taq (Tli RNaseH 
Plus, Takara Biotechnology, Japan), 0.4 uL Forward Primer, 0.4 uL Reverse Primer, 
0.4 uL Rox Reference Dye II (Takara Biotechnology, Japan), 2 uL DNA template and 
6.8 uL ddH20. Standard curves were obtained by serially diluting (10? - 10? copies 
ИГ!) already quantified plasmid DNA containing the AOB-amoA, AOA-amoA, nirK or 
nosZ gene fragment. The standard plasmids were prepared from the extracted DNA 
samples using same primers as above. The PCR amplification efficiency - R? was 
103.0% - 0.997 for AOB-amoA, 101.0% - 0.982 for AOA-amoA, 99.8% - 0.980 for 
nirK and 95.3% - 0.994 for nosZ. Gene reaction programs were listed in Table S2 


based on previous studies (Rich et al., 2003; Levy-Booth et al., 2014). 
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2.4 Statistical analyses 


One-way ANOVA analyses followed by least significant difference (LSD) tests 
were used to test the significant differences of soil properties, N2O emission rates, 
gross N transformation rates and functional gene abundances among different N 
addition levels at each site. Two-way ANOVA analyses were performed to compare 
the difference of N2O fluxes, N transformation rates and gene abundances between the 
short-term and long-term sites. Two pairs of N addition levels at the short-term site 
(SL: 35 kg М ha! а"; SH: 105 kg М па a!) and the long-term site (LL: 50 kg М ha’! 
а; LM: 100 kg М ha'! а!) were compared correspondingly. Before one-way and two- 
way ANOVAs, data normality was tested and logarithmic transformation was applied 
if necessary. To quantitatively compare the N addition effects between the two sites, 
we also calculated the relative change percentage (%) of soil N2O fluxes, N 
transformation rates and gene abundances under N addition treatments to those under 
the control treatment. Moreover, Pearson correlation analyses (SPSS 17.0, SPSS Inc., 
Chicago, USA) were conducted to detect the relations among soil properties, gene 
abundances and №0 fluxes for both sites. Statistically significant differences were 
analyzed for P < 0.05 values unless otherwise stated. All statistical analyses were 


performed in SPSS 17.0 (SPSS Inc., Chicago, USA). 


3. Results 


3.1 Background soil chemical properties at the short- and long-term sites 


Soils from both sites were strongly acidic with average pH of 3.8 at the long-term 
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site and that of 3.9 at Ше short-term site (р < 0.05, Table $3). Soil NO3-N 
concentration was significantly higher than NH4'-N concentration in the soils for both 
sites, while the short-term site had a relatively lower NO3-N content compared to the 
long-term site. The long-term site soils also had relatively higher TOC content (p « 
0.05) and higher total-N content than the short-term site soils (Table S3). No 
significant difference for the C/N ratio was found among the four N treatments within 


each site and between the two sites (p > 0.05, Table S3). 


3.2 N concentrations and УМ enrichments in inorganic-N and МО 


During the 14 days of tracing incubation, soil NH4* concentration decreased on 
average by 7.7% and by 15.5% across the four М treatment levels for the short- and 
long-term sites, respectively; whereas soil МОз` concentration increased on average 
by 25.0% and 10.2% for the short- and long-term sites, respectively (Fig. S2a, b, d, e). 
A decline in the ^N enrichment of the NH4'-N pool was observed when ammonium 
was labelled (Fig. S2c), indicating that natural abundance or low abundance of NH4* 
entered this pool. A similar declining trend was observed in Ше РМ enrichment of the 
NO-N pool when nitrate was labelled (Fig.S2f), showing natural ог low abundance 
of NO; entered the pool. Furthermore, the ^N enrichment of the МОз-М pool 
increased in the ammonium labelled treatment (Fig. S2c), indicating a conversion of 
labelled МН“ to NOs (e.g. ammonia oxidation). The modelled and observed data 
agreed well with the R? being generally larger than 0.92. 

The '°N enrichment in №0 ranged from 0.01 to 0.53 atom% excess during the 
incubation, and showed initially an increase followed by a decreasing trend towards 


the end of the incubation (14 d) (Fig. Та, c). In addition, №0 emitting from the nitrate 
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labelled soils showed а significantly higher ^N enrichment than N20 emitting from 
the ammonium labelled soils during the incubation (р < 0.05, Fig. Та, с). A 
significantly higher УМ О/РМОУ % was found compared to the ^N;0/?NH4'96 


throughout the whole incubation (Fig. 1b, d). 


3.3 №0 fluxes and production pathways 


With increase in the level of М addition, soil average №0 fluxes of the low or 
medium М treatments during the incubation also increased, while N2O fluxes of the 
high М treatments decreased at the two sites (Fig. 2a). The N2O flux in SL, SM and 
SH was 42.3%, 2.7% higher and 11.6% lower than the treatment control (SC) at the 
short-term site, while the ХО flux in LL, LM and LH was 35.7%, 378.5% and 
106.7% higher than the control (LC) at the long-term site. Except for the significantly 
higher flux in LM, №0 fluxes were not significantly different among М treatments or 
between sites (p > 0.05). As to №0 emission pathways, denitrification (d) contributed 
more than 40% to soil №0 production, while the trend of this pathway within 
increasing М levels differed at the two sites (short-term: SC<SH<SL<SM; long-term: 
LL<LM<LC<LH) (Fig. 2b). The contribution of heterotrophic (h) and autotrophic 
nitrification (a) was in a range of 17.7% - 36.8% and 14.3 - 23.2%, respectively. 
However, the two nitrification fractions of N2O differed at the short- vs. long-term 
within increasing N levels (h, short-term: SM<SC<SL<SH, long-term: 
LC<LH<LM<LL,; a, short-term: SL<SH<SM<SC, long-term: LM<LL<LH<LC) (Fig. 
2b). Furthermore, the contribution of different N2O emission fractions changed from a 
dominant heterotrophic nitrification (N2Ox) to denitrification (№Оа) towards the end 


of the study, while the contribution of autotrophic nitrification (N20a) increased in the 
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early stage (0 - 7d) of the incubation (Fig. 2c). The contribution of the three fractions 
was on average 53% (d), 28% (h) and 18% (a) during the incubation, showing an 


overall dominance of denitrification in the studied subtropical forest soils (Fig. 2c). 


3.4 Gross and net N transformation rates 


Gross autotrophic nitrification rate (Омна) was highest among all the ten М 
transformation rates under almost all N treatments, followed by gross mineralization 
rates of recalcitrant and labile organic N to NH4* (Муес and Ммађ) (Table S4). Gross 
heterotrophic nitrification rate (Омес) was much lower than gross autotrophic 
nitrification rate in the studied subtropical forest soils (Table S4). Total gross 
nitrification rate (Охна + Омес, not shown in the table) ranged from 0.08 mg kg"! d`! to 
0.19 mg kg! d! while Омнч contributed over 99% to the total NOs' production. 
Except for SC and LH treatments, soil gross NO3 immobilization rate (Гуоз) was 
higher than gross NH4* immobilization rates (Гунч-меез, IvH4-Nlab) and gross rates of 
dissimilatory NOs reduction to NH4* (ОМКА) (Dyo3) under most N treatment 
conditions (Table S4), indicating that МОз was more easily trapped in the organic М 
pool than in ће NH4'-N pool. Net NOs' production rate (МЕТуозрюа) was much 
higher than net NH4* production rate (VETwo3proa.), while LH had the highest net 
МНа loss but the lowest net NOs production from all the eight М treatments. 

МЕ Туозргод. also differed significantly among the four N-treatments in either the short- 
term or the long-term site (p < 0.05). АП the ten gross М transformation rates 
considered in both sites were statistically different with respect to the N treatments (p 
< 0.05, Table S4). However, according to two-way ANOVA results, long-term М 


addition did not significantly change М transformation rates in soils (p > 0.05) in 


18 


comparison to Ше short-term site. 

The response direction of gross organic М mineralization rates (Ммаь and Mnrec) 
to N addition was opposite between short-term and long-term sites (Fig. 3a, b), so did 
the gross rate of heterotrophic nitrification (Ол) (Fig. Зе). Gross mineralization rate 
of labile organic N to NH4* (Ммаь) was 82%, 144% and 58% higher respectively in 
SL, SM and SH than that in SC at the short-term site, but the Мма was 44%, 20% and 
47% lower respectively in LL, LM and LH than that in LC at the long-term site (Fig. 
За). In contrast, the gross mineralization rate of recalcitrant organic N to NHa4* (Mnrec) 
was 49%, 16% and 15% lower respectively in SL, SM and SH than that in SC at the 
short-term site, whereas the Мугс was 51%, 85% and 128% higher respectively in LL, 
LM and LH than that in LC at the long-term site (Fig. 3b). The large stimulation 
effect of М addition on the gross rates of NO3 immobilization (Inos) and dissimilatory 
NOs reduction to NH4* (ОМКА) (Dwos) at the short-term site, decreased or even 
became inhibition at the long-term site (Fig. 3d, В. Compared with the short-term site, 
the high N addition treatment at the long-term site induced the largest inhibition effect 
on the gross autotrophic nitrification rate (Омна) and the net МОз production rate 


(NETwo3-prod) (Fig. 3c, g). 


3.5 Functional gene abundance 


AOA-amoA was the most abundant, followed by nirK, nosZ and AOB-amoA, 
while АОВ-атоА was least abundant at both the short-term and long-term sites (Fig.4 
а, b, d, e). The ratio between denitrifying and nitrifying gene abundance, i.e. (nirK 
+nosZ)/(AOB-amoA + AOA-amoA) was mostly less than 1.0 at both sites (Fig. 4c), 


indicating that nitrifiers were more abundant than denitrifiers in the forest soil. The 
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nosZ/nirK ratio was mostly less than 1.0 а both sites (Fig. 4f), indicating Ша! N20- 
producers were less abundant than N2O-reducers. Although AOA-amoA, nirK and 
(тїї К +nosZ)/(AOB-amoA + AOA-amoA) ratio responded significantly to N addition 
within each site (p < 0.01) (Fig. 4b, c, d), the abundance of all assessed functional 
genes showed no significant difference between the two sites (p > 0.05). 

Similar stimulation effects of М addition on nirK abundance were observed at the 
two N addition sites (Fig. 4d, inset). Except the high N treatments of the two sites, the 
direction of nosZ abundance, (nirK +nosZ)/(AOB-amoA + АОА-атой) and 
nosZ/nirK ratios was also the same (Fig. 4 с, е, Е, insets), indicting a stimulation effect 
caused by low and medium N additions in the forest soil. However, except the high N 
treatments of the two sites, the direction of AOB-amoA abundance response to М 
addition was opposite between the two sites (Fig. 4a, inset). The abundance of AOA- 
атоА was 26.5% lower, 72.0% higher and 127.3% higher respectively in SL, SM and 
SH than that in SC at the short-term site, while the AOA-amoA abundance was 46.8% 
lower, 8.0% higher and 27.7% lower respectively in LL, LM and LH than that in LC 
at the short-term site (Fig. 4b, inset), indicating that under high N treatment the AOA- 
amoA abundance was stimulated at the short-term site but suppressed at the long-term 


site. 


3.6 Correlations among soil properties, abundance of functional genes and N20 


fluxes 


The correlation relations among soil properties, functional gene abundance and 
№0 fluxes differed between short-term and long-term sites (Fig. 5). At the short- 


term site, soil №0 flux was positively correlated with nosZ abundance (p < 0.05). The 
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abundance of nosZ was also positively correlated with both nirK (р < 0.05) and AOB- 
атоА abundances (p < 0.01), implying there might exist possible interactions between 
nitrifying and denitrifying groups in regulating N20 emission. AOA-amoA abundance 
was negatively correlated with soil pH (p « 0.01) but positively with soil C/N ratio (p 
« 0.05) (Fig. 5). The abundance of AOB-amoA and nosZ was positively correlated 
with soil NH4'-N, total N and C/N ratio (p < 0.05) (Fig. 5), indicating that nitrifying 
N2O-producers and denitrifying N2O-reducers were also linked to substrate 
availability in the acid soil environment. Positive correlations were also found 
between (nirK +nosZ)/(AOB-amoA + AOA-amoA) ratio and soil NH4'-N and 
between nosZ/nirK ratio and TOC (Fig. S3a), further indicating that the composition 
of nitrifiers and denitrifiers could possibly be affected by soil substrate availability. 
However, at the long-term N addition site, positive correlations were only found 
between №0 flux and AOA-amoA abundance (p < 0.05), and between AOB-amoA 
abundance and pH (p « 0.05) (Fig. 5). In general, the correlations among soil 
properties, functional gene abundance and N20 fluxes was tighter at the short- than 


the long-term site. 


4. Discussion 


4.1 Responses of N20 fluxes and МО production pathways to short- and long-term N 


additions 


The stimulation effect of N addition on ХО emission rates was obvious only in 
the low and medium М treatments at the two sites (Fig. 2а). The №0 production- 


consumption balance (Xu et al., 2012; Zhang et al., 2014) during the denitrification 
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process of the forest soil may be changed: low and medium М additions favored N20 
production but not №0, whereas this effect became weaker under high М addition. 
The more dominant №0 production than consumption is also reflected by the more 
abundant N2O-producer (nirK) than N2O-reducer (nosZ) (Fig. 4 d, e, f), which implied 
a higher activity of nirK than nosZ in the subtropical forest soil (Zhang et al., 2014). 
However in the high М treatments, N20 emissions were reduced at the short-term site 
but still stimulated at the long-term site compared to the controls (Fig. 2a). This is in 
contrast to the results from a montane forest with 3-4 years of N addition compared to 
11-12 years N addition of a lowland forest (Corre et al., 2014). Possible reasons might 
be that soil microbes responsible for №0 emissions were more sensitive in 
responding to short-term and higher N additions, but gradually adapted over longer 
time periods, so that the stimulation effect at high-N subsided. Although long-term N 
addition did not induce significant changes in total N2O emission, the contribution of 
the various emission pathways changed, with higher denitrification contribution 
(53%) and lower heterotrophic (28%) and autotrophic nitrification contributions 
(18%) (Fig. 2b). Thus denitrification remained the dominant pathway for N2O 
emissions in this subtropical forest soil, which is consistent with the conclusion of 
Zhang et al. (2011). The ^N;0/^NH4'96 and !°N20/!°NO3% were not significantly 
different among increasing N levels or between short-term and long-term sites (Fig. 
Ib, d). This is not in line with previous findings which showed that increasing 
ammonium concentrations supported a higher contribution of nitrification to N20 
production (Müller et al., 1998; Avrahami et al., 2002). Possible reasons might be that 
the top soils from this subtropical forest had been naturally N-saturated in the control 
treatments (SC and LC) (Gurmesa et al., 2016), thus causing a non-recognizable 


response of the nitrification or denitrification contribution to additional NH4* or ХО». 
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By combining our observed М rate dynamics, we suspect that the added NH4* in Ше 
forest soil was quickly transformed into МОз, and subsequently denitrified into 
gaseous М (N20). This coupled nitrification-denitrification process is prevalent when 
rates of autotrophic nitrification are high and may have contributed to the low 
contribution of nitrification related N2O (Müller et al., 2014). This process may also 
contributed to the significantly lower РМ enrichment in N20 in contrast to the higher 
PN enrichments іп МНА? and NO: in the subtropical forest soil (Fig. Та, c; Fig. S2 c, 
f). 

We separately labelled the NH4*-N pool and NOs-N pool in the paired ^N 
tracing incubation. The ^N enrichment in NO3 was higher than the PN enrichment in 
№0 in the РМНА ЧађеПеа treatment. The PN enrichment in NH4* was also higher 
than the ^N enrichment in N20 in the УМОз labelled treatment (Fig. S2). This 
further implies that the transformation from NH4” to NOs' was more prevalent than 
the transformation from inorganic N to ХО. Autotrophic nitrification only 
contributed 18% to the total №20 flux, despite the dominant status of this process in 
the ten N transformation processes which supports the above reasoning. This may 
indicate that fast autotrophic nitrification induced №0 production may also have 
stimulated further №0 reduction to №, possibly via nitrifier-denitrification (Müller et 
al., 2014). Thus, the measurement of № and its PN enrichment would further 
elucidate this possible process. Heterotrophic nitrification contributed considerably 
(28%) to №0 fluxes while only a tiny part (1%) contributed to total NO3 production 
(Table S4, Fig. 2c). This could be explained by the “hole-in-the-pipe” model and N20 
gas leakage characteristics in subtropical forest soils: the higher autotrophic 
nitrification rate induced a bigger “hole” for №0 leakage while the lower 


heterotrophic nitrification rate induced a smaller “hole” for №0 leakage, thus 
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inducing Ше comparability of Ше two №0 production pathway contributions 


(Firestone and Davidson, 1989; Zhang et al., 2014). 


4.2 Responses of soil N transformation rates to short- and long-term N additions 


Except for the high N treatment at the long-term site, the subtropical forest soil 
which received both natural (SC, LC) and anthropogenic N depositions showed higher 
gross autotrophic nitrification and mineralization rates than other N transformation 
rates (Table S4). According to Cheng et al. (2014), the relative fast transformations 
from КН“ to NO3 and from organic М to NH4* might be one of the reasons for the 
natural N enrichment in the humid subtropical forest soil. Moreover, the more 
dominant autotrophic nitrification than mineralization processes (except for the high 
М treatment at the long-term site) also well explains the higher NO3 than NH4* 
concentrations in the forest soil (Table S3). The particularly higher gross МНА? 
immobilization rates (Гуна-мгес, [мн4-Моь), higher gross mineralization rates (Мус) but 
lower gross autotrophic nitrification rates (Own) in the long-term high М treatment 
(LH) (Table S3) suggest that the soil mineralization-immobilization turnover (MIT) 
became more dominant than NH4*-NO3° turnover (i.e. autotrophic nitrification) after 
long-term anthropogenic N addition. In addition, the significant differences of soil N 
transformation rates among N treatments at both sites (Table S4) further strengthens 
the implication importance of N transformation responses for the N enrichment, N 
loss and soil acidification in subtropical forest ecosystems (Gao et al., 2015). 
However, N transformation rates were not different significantly between the short- 
term and long-term sites (Table 54). Previous results revealed that soil М 


transformation rates are affected by a series of factors, including the interacted N 
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transformation processes (nitrification-mineralization, mineralization- 
immobilization), soil microbes (microbial activity and functional abundances) and soil 
properties (C/N ratio, TOC and pH) (Bengtsson et al., 2003; Zhu et al., 2013b; Gao et 
al., 2015; Chen et al., 2017). So the insignificant change of soil N transformation rates 
at the long-term site might be owed to that, although long-term anthropogenic N 
addition still induced significant changes in soil N transformation rates, the interacted 
soil microbes which play roles in N transformation processes have adapted to current 
soil conditions. The combined effects of diverse soil abiotic and biotic factors on N 
transformations were weakened, thus finally preventing the possibly large responses 


of N transformations to prolonged N addition. 


Our results showed that the rate responses of specific N transformation processes, 
which mainly associate to the fate of NH4* and NOs, were different between the 
short- and long-term sites (Fig. 3). Previous studies have revealed that soil organic 
matter fractions are linked to different functional microbes while their responses to N 
depositions are also different (Frey et al., 2004; Cusack et al., 2011). So the 
quantification of the gross mineralization rates of labile organic М (Ммаь) and 
recalcitrant organic М (Мег) will facilitate the understanding of the fast (labile 
organic N) and slow (recalcitrant organic N) turnover of soil organic N pools, 
respectively (Gao et al., 2016). In our study, Ммаь was stimulated by short-term М 
addition but was inhibited by long-term N addition, while results in Му are totally 
opposite (Fig. 3a, b insets), showing different response directions of gross 
mineralization rates not only between labile and recalcitrant organic N forms but also 
between short- and long-term N addition treatments. Possible reasons might be that, 
short-term N addition limited the biomass and activity of the soil fungi which 


participant into the mineralization of recalcitrant organic N, thus indirectly inducing 
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an inhibition in Мус (Frey et al., 2004; Schimel and Bennett, 2004). Moreover, М 
addition may also enhance the transformation of soil labile organic matter to 
recalcitrant organic matter (Maaroufi et al., 2015). This might be one of the reasons 
for the stimulated Ммав but the inhibited Ммав by short-term М addition. The 
stimulation effect of long-term N additions on Myrec might be owed to that, the 
microbial N limitation was alleviated after long-term N addition, thus inducing an 
increased microbial activity for the stimulation of organic matter decomposition 
(Koranda et al., 2014; Gao et al., 2016). The inhibited Mya at the long-term site 
might be a downregulation response to the microbial upregulation to Му, while the 
opposite responses in Ммиес and Ммар under М additions could be contributed by the 
changed stabilization of soil organic matter into recalcitrant compounds, and also by 
the changed microbial extracellular enzyme activity which conduct the production and 
decomposition of soil organic compounds (Sinsabaugh et al., 2005; Janssens et al., 
2010; Gao et al., 2016). The response of soil gross autotrophic nitrification rate was 
most significant, except for the high N treatment of the long-term site (Fig. 3c inset), 
indicating that under this highest М addition rate (i.e. 150 kg М ha! a’), soil may have 
changed most and thus contributed much less to nitrifying №0 emission. The special 
change of soil in the high N treatments of the long-term site could also be reflected by 
the lowest net NO3 production rate and the highest NH4* immobilization rate in 
comparison to other N treatments (Table S4). Soil gross heterotrophic nitrification 
rate (Onrec) was stimulated by short-term М addition but was inhibited by long-term М 
addition (Fig. За), implying that the production of №0 and NOs via heterotrophic 
nitrification pathways tented to increase at the short-term but decrease at the long- 
term site. The more positive responses of Омиес, [моз and Dyos at the short-term than 


long-term site prevent a higher loss risk of soil NO3 through gaseous М emission 
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(denitrification) (Wang et al., 2015). 


4.3 Regulations of soil nitrifiers and denitrifiers in N2O emissions and their abiotic 


controllers 


Previous study exhibited higher abundance of ammonia-oxidizing archaea (AOA) 
in the same subtropical forest which is consistent with our present finding (Isobe et 
al., 2012). Moreover, several studies have highlighted a predominant role of AOA in 
nitrification (Gubry-Rangin et al., 2010; Verhamme et al., 2011; Zhang et al., 2012; 
Faeflen et al., 2016). So the dominant autotrophic nitrification in our current results 
(Table 54) might be explained by the more abundant AOA than AOB in the strong 
acidic soil (Zhang et al., 2012; Faeflen et al., 2016). For a better defining the relative 
importance of AOA and AOB in the nitrification processes, future explorations of the 
archaeal and bacterial ammonia oxidizer communities needs to be considered (Gubry- 
Rangin et al., 2010). Functions of nirK and nosZ in relationship to denitrification 
processes have proved to be different: nirK plays roles in N2O production while nosZ 
play roles in №0 reduction or consumption (Zhang et al., 2014; Domeignoz-Horta et 
al., 2017). Our results revealed a slightly higher abundance of nirK than nosZ (Fig. 
Ad, e, f), potentially indicting that №0 production exceeds N20 reduction and further 
contributing to N2O emission out of the forest soil. This is also in line with previous 
findings that the studied subtropical forest soil plays as a natural source of №0 (Fang 
et al., 2008а; Tang et al., 2006). However, it is worth mentioning that NO as another 
intermediate of denitrification cannot be neglected in considering gene regulation 
functions since nirK plays essential roles in the second step of denitrification, i.e. 


nitrite reduction to NO (Levy-Booth et al., 2014). Moreover, nirS as another 
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molecular marker gene of denitrifying bacteria for this step still needs to be fully 
understood (Prieme et al., 2002; Barta et al., 2010). The nosZ/nirK ratio in most М 
treatments was less than 1 (Fig. 4f) and finally led to an average value of 0.68 (data 
not shown) without considering the effects of N addition. This is dramatically 
opposite with the findings of Zhang et al. (2014) which uncovered an average value of 
0.7 for nirK/nosZ ratio in subtropical soils. Possible reasons might be that, the strong 
acidic as well as highly weathered soil in our subtropical forest contributes to a higher 
redox potential, which supports a more active nirK than nosZ communities in the soil 
(Qafoku et al., 2004; Xu and Cai, 2007; Zhang et al., 2014). 

Compared with the short-term N site, the prolonged N addition at the long-term 
site did not induce significant changes in functional gene abundances, but changed the 
responses of specific functional gene abundances to N addition treatments markedly 
(Fig. 4). This further highlights a higher importance of functional gene response 
patterns to N addition levels than gene responses to N addition durations. Possible 
reasons for the insignificant changes of gene abundances under the long-term N 
additions might be that, the activity of soil nitrifiers and denitrifiers are affected by 
soil abiotic conditions to a large extent, including pH, redox potential, N availability 
and other nutrients (Avrahami et al., 2002; Barta et al., 2010; Hu et al., 2013; Yuet al., 
2014; Zhang et al., 2014). Thus the comprehensive effects of soil abiotic factors 
eliminated possibly significant responses of functional groups with the N addition 
duration prolong. We observed significant responses of AOA-amoA and nirK 
abundances but not AOB-amoA or nosZ abundances to М addition treatments at the 
two sites (Fig. 4). This suggests different response sensitivities of soil functional 
microbial communities to N addition (Avrahami et al., 2002; Li and Gu, 2013). 


Significant positive correlations were observed between AOB-amoA and nosZ 
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abundances, and also between nirK and nosZ abundances at Ше short-term site (Fig. 
5), indicating a combined gene regulation on №0 emissions from both nitrification 
and denitrification pathways (Zhang et al., 2014; Tang et al., 2016; Domeignoz-Horta 
et al., 2017). However, this close relationship linking both nitrifying and denitrifying 
functional groups was not found at the long-term site. Possibly because that the short- 
term N addition induced more significant changes in soil properties but this effect was 
weakened by the long-term N addition. This could be explained by the more 
significant correlations between soil properties (pH, C/N ratio, total М and NH4'-N) 
and functional gene abundances at the short- than long-term site (Fig. 5). The positive 
correlation between nosZ abundance and №0 flux observed at the short-term М 
addition site seems a contradiction, possibly because the N2O flux we calculated was 
a sum of three N20 emission fractions associating autotrophic nitrification, 
heterotrophic nitrification and denitrification pathways (Zhang et al., 2011). However, 
at the long-term site, we only found a significant negative correlation between pH and 
AOB-amoA abundance, and a significant positive correlation between AOA-amoA 
abundance and N20 flux (Fig. 5), implying a relatively looser relationship among soil 
properties, gene abundance and №0 fluxes. Previous study suggested that soil 
properties and microbial community may distinctly play dominant roles in explaining 
lower and higher N20 emission rates, respectively (Domeignoz-Horta et al., 2017). 
This triggers our speculation of possible reasons for the quite different relationships of 
soil, gene and N20 between the short- and long-term sites: the deepened acidification 
in the long-term site soil induced not only the community changes but also the 
ecological niche distinctions of the two ammonia oxidizers (AOA and АОВ) (Nicol et 
al., 2008; Hu et al., 2013; Remy et al., 2018). Moreover, the change of soil conditions 


under the long-term N addition also reduced the close relationship between AOA and 
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other variables in this acid soil (Hu et al., 2013) and finally induced Ше positive 
correlation between pH and AOB-amoA abundance. Briefly, the changed soil and 
climatic conditions (rainfall, moisture, temperature, etc.) affected microbial activities 
and gene abundance related to N2O production and further reduction to N2 which is 
affected by soil acidity (Liu et al., 2014). Thus, the relationship among soil properties, 
gene abundance and N20 fluxes at the short-term N addition site is much tighter than 
that at the long-term site where microbial communities have adapted to prevailing 


conditions. 


Our results showed a changed functional gene regulation in the subtropical forest 
N2O emission as a response to prolonged N addition duration. However, it is still 
obscured for the understanding of soil nitrifying and denitrifying functional group 
changes from only the functional gene abundance aspect since the operating of gene 
function is reflected by corresponding gene transcripts (Xu et al., 2012; Wu et а1., 
2017). So for a better elucidating the importance of the internal mechanisms of soil 
microbial functional groups in regulating N2O emissions, studies of corresponding 
RNA transcripts and their relationships with N2O productions (Liu et al., 2014; Liu et 
al., 2010; Hassan et al., 2016) are needed to be conducted in the future. Moreover, the 
comprehensive exploration of soil microbial responses to N depositions from 
microbial structures, gene abundance and transcription levels need to be done in the 


future (Li and Gu, 2013). 


4.4 №0 productions in relation to М additions, N transformations and functional gene 


regulations 


Consistent with previous results (Zhang et al., 2011), the contribution of 
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heterotrophic nitrification was lower than denitrification but was higher than 
autotrophic nitrification to total N2O fluxes (Fig. 2b). However, the higher 
heterotrophic than autotrophic nitrification contribution is not reflected by either the 
significantly higher gross autotrophic nitrification than heterotrophic nitrification rates 
or by the most abundant AOA-amoA (Table S4, Fig. 4b), while the most dominant 
denitrification contribution is also inconsistent with the relatively less abundant 
denitrifiers (nirK + nosZ) than nitrifiers (AOA-amoA + AOB-amo4A) in the forest soil 
(Fig. 4c). Based on the conceptual “hole-in-the-pipe” (HIP) model of N20 emission 
and the modifications of the model (Firestone and Davidson, 1989; Zhang et al., 2014; 
Zhang et al., 2015), we speculate possible reasons for the inconsistence among 
functional gene abundance, N transformation rates and N2O emissions: the more 
abundant AOA in soil contributed to a higher autotrophic nitrification rate but also 
triggered a bigger N2O leaking “hole” from the “pipe” of this pathway. The lower 
heterotrophic nitrification rate induced a smaller “hole” for №0 leaking, so a 
considerable proportion of N2O production via this pathway maintained. 
Denitrification plays as an essential pathway for NO3 removing from this subtropical 
forest (Fang et al., 2015), but the higher nirK than nosZ abundance in the forest soil 
potentially induces more N20 production than №0 reduction, thus finally induced the 
highest contribution of denitrification to N2O emission (Fig. 2b). 

Our study in the subtropical forest revealed a change of the response patterns in 
specific N transformation rates to long-term N addition (Fig. 3), further directing a 
change in the transformation fates of NH4* and МОз, including the production of 
МН via the mineralization of soil labile and recalcitrant organic N and the 
consumption of NH4” via autotrophic nitrification, the production of NOx via 


autotrophic and heterotrophic nitrifications and the consumption of NO3 via 
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immobilization апа ОМКА, respectively, to Ше prolonged М addition. Thus Ше 
different responses of net МОз production to increased N addition at the two sites 
(Fig. 3g) further reflect a change of NO3 retention potential in the subtropical forest 
soil with N addition duration prolong. Our study revealed that soil functional genes, 
especially the two nitrifying functional genes (AOB-amoA, AOA-amo4A) changed 
more in responding to prolonged N addition (Fig. 4). However, it is still controversial 
to conclude that the insignificant change of N2O under the long-term N addition is a 
combined response pattern of the four genes that currently studied. The emission rate 
of №0 is a reflection balancing the regulations among nitrifying N2O-producers 
(AOB-amoA, AOA-amo4A), denitrifying N2O-producers (nirK/nirS) and denitrifying- 
№0 reducers (nosZ). So the fully understanding of other functional groups that not 
considered in the present study but are of great importance for N2O emission, 
including nirS gene for №0 production from NO, and the newly identified N20- 
reducting clade recently for the N2O reduction to N2, needs to be done in the future 
(Levy-Booth et al., 2014; Zhang et al., 2014; Tang et al., 2016; Domeignoz-Horta et 
al., 2017). Correspondingly, NO as an intermediate production and Хо as a final 
production in both nitrification and denitrification processes cannot be neglected in 
future studies because of their close linkage with N2O production and consumption 
(Zhang et al., 2014; Zhang et al., 2015; Tang et al., 2016). 

In this study, we simulated М deposition by adding equal amounts of NH4* and 
NO (МНАМОз) to subtropical forest soils. Natural atmospheric М depositions usually 
show unequal contributions from these two inorganic N forms, with a ratio of 
approximately 60/40 (NH4'/NOs) (Fang et al., 2011; Fang et al., 2015). A study of 
two subtropical masson pine forests in south China even found a ratio of about 70/30 


(Chen et al., 2004), thus providing different transformation capacities of deposited 
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МН and NOs. The studied subtropical forest in DHSBR has received atmospheric N 
deposition at an amount of more than 30 kg М ha’! a`! since 1990, so the long-term 
natural N accumulation in the forest soils might be one of the reasons for the 
significant net NO3 production in our findings, although leaching and denitrification 
have been proved to be main fates of the produced МОз` in this forest (Fang et al., 
2009; Fang et al., 2015). Although high autotrophic nitrification rates were measured 
in almost all the N addition treatments at the two sites, we did not find a 
correspondingly high N enrichment іп №0. One reason for this result could be high 
conversion rates of N20 to №. So the determination of the ^N enrichment in № 
should be considered in future studies. Despite the obvious changes of specific soil N 
transformation rates and gene abundances but not N2O fluxes in responding to long- 
term N addition, we still found a looser relationship among soil properties, functional 
gene abundance and N20 fluxes at the long-term site than the short-term site based on 
correlation analysis results. We argue that functional microbial community 
adaptations occur which are developing in response to prevailing soil conditions and 
high natural ambient N depositions. This is in line with the response of microbial 
community structure, diversity and also gene expression to changing biotic and 
abiotic factors (Avrahami et al., 2002; Smith et al., 2010; Németh et al., 2014), while 
the effect of N deposition is a prevailing factor over other influencing factors. This 
would also explain the change of the correlations between functional gene abundance 
and №0 flux with the М addition duration prolong (Fig. 5). However, to elucidate this 
further microbial activity shifts rather than just focusing on gene expressions would 


most likely be more indicative. 


5. Conclusion 
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In this study, we compared Ше response of №0 emissions, N transformations as 
well as gene regulation traits between two subtropical evergreen broadleaf forest sites 
that have received either one year (short-term) or 12 years (long-term) of simulated N 
deposition. Our results revealed that, long-term N addition did not induce significant 
changes in №0 fluxes, soil М transformation rates and functional gene abundances. 
However, compared with the short-term N addition: i) long-term N addition changed 
the response of specific N transformation rates, especially the rate of N transformation 
processes that related to the fate of МНА and NOs (mineralization, nitrification, ХОз 
immobilization and the dissimilatory NO3 reduction to МН4*); ii) long-term М 
addition induced more changes in the response of soil nitrifying N2O-producters 
(AOA-amoA, AOB-amoA) to М treatments than that in the response of denitrifying 
N2O-producers (nirK) or N2O-reducers (nosZ); ш) long-term N addition weakened Ше 
close correlations among soil properties, functional gene abundance and №0 fluxes 
that observed at the short-term site. However, the observed differences of N dynamics 
and №0 emission traits at the two simulated М deposition sites cannot be solely 
explained by soil microbial nitrifying and denitrifying functional genes. Soil 
microenvironment conditions (pH, substrate availability, etc.) and their transition 
effects from the short-term to long-term N depositions are key to understanding 
microbial functional gene expressions and gene associated N dynamics (including 
№0). We argue that microbial regulation adaptations exist to prevailing soil 


conditions in response to long-term natural and anthropogenic N depositions. 
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Figure legends 


Fig. 1. The PN enrichment in №0 during the paired PN tracing incubation (а, с), and 
ООЉРМНА ratios in the ammonium-labelled incubation (b) and ^N20/PNOX* ratios 
in the nitrate-labelled incubation (d). Symbols (points and lines) in light green and 
orange denote the short-term N addition treatments and long-term N addition 
treatments, respectively. SC, SL, SM and SH denote control, low-N, medium-N and 
high-N treatments at the short-term N addition site, respectively. LC, LL, LM and LH 
denote control, low-N, medium-N and high-N treatments at the long-term N addition 
site, respectively. *, **, ***: the difference between PNH4'-labelled treatment and 
I5NO;--Iabelled treatment is significant at p < 0.05, p < 0.01 and p < 0.001 level, 


respectively. 


Fig. 2. Soil №0 average fluxes (a), specific pathway fractions of emitted N2O in each 


N treatment (b) and average pathway contributions over the 14 day incubation of the 
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subtropical forest soils (с). SC, SL, SM and SH denote control, low-N, medium-N and 
high-N treatments at the short-term N addition site, respectively. LC, LL, LM and LH 
denote control, low-N, medium-N and high-N treatments at the long-term N addition 


ВЕЗИ 


site, respectively. The sign above the light green bars in panel a denotes the 
difference between the corresponding N treatment and other N treatments is 
significant (p < 0.05). Bars in light green, moderate green and dark green in panel b 
denote the contribution of N2O fraction that originate from heterotrophic nitrification 
(h), autotrophic nitrification (a) and denitrification (d), respectively. Areas in light 
gray, moderate gray and dark gray in panel c denote that during the incubation, the 


contribution of N20 fraction that originate from heterotrophic nitrification (NO), 


autotrophic nitrification (N2O;) and denitrification (N2Oq), respectively. 


Fig. 3. The PN tracing model (Müller et al., 2007) and relative change percentages of 
soil gross/net N transformation rates in the N addition treatments in comparison to 
controls at the two sites (inset a-g). White squares and black arrows indicate soil N 
pools and М transformation processes in the tracing model, respectively. NHa4*: 
ammonium; NHa*aas : adsorbed МНА”; Nia»: labile organic N; NOs’: nitrate; Мес: 
recalcitrant organic М. M rec: mineralization of recalcitrant organic nitrogen (о МНА”; 
Inna-Nrec: immobilization of NH4* to recalcitrant organic nitrogen; Ммаь: 
mineralization of labile organic nitrogen; Гуна-мађ: immobilization of МНД? to labile 
organic nitrogen; Onrec: oxidation of recalcitrant organic nitrogen to МОз ; Ios: 
immobilization of NO: to recalcitrant organic №; Onpa: oxidation of NH4* to МОз ; 
Dyos: dissimilatory NO3 reduction to NH4*; Ayna: adsorption of NH4* on cation 
exchange sites; Куна: release of МН“ on cation exchange sites; NETwosaproa: net 
NO» production. Insets numbered with a, b, c, d, e, f and g denote the relative change 
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percentage of Ммаь, Мугес, Омна, Омоз, Омгес, [моз, and МЕТмозарта, respectively. 
Within each inset, bars in light green and orange denote Ше results at Ше short- and 
long-term sites, respectively. L, M and H in the x-axis denote the low-N, medium-N 
and high-N treatment, respectively. Numbers in the y-axis denote the relative change 
percentage (%) of individual N transformation rate in the N addition treatment in 


comparison to control. 


Fig. 4. Abundance (copy numbers g`! dry soil) and abundance ratios of the functional 
genes (AOB-amoA, АОА-атоА, nirK and nosZ) under different М treatments at the 
two sites. Bars in light green and orange denote short-term N addition treatments and 
long-term N addition treatments, respectively. Inset denote relative change percentage 
of gene abundance or abundance ratio in N-addition treatments in comparison to 
control. L, M and H in the x-axis denote the low-N, medium-N and high-N treatment, 
respectively. Numbers in the y-axis denote the relative change percentage (%) of gene 
abundance or abundance ratio in N addition treatments in comparison to control. D 
and N indicate results of two-way ANOVA analysis: D indicates the difference 
significance between М addition durations (short-term and long-term), and N (М level) 
indicates the difference significance among N treatments. Only significant differences 


were showed in this figure while insignificant differences were not shown. 


Fig. 5. Pearson correlation analysis results of soil properties, gene abundances and 
№0 fluxes at the short-term and long-term sites. N2O production pathways that 
considered in this study, i.e. heterotrophic nitrification, autotrophic nitrification and 
denitrification, and the roles of the studied functional genes (AOA-amoA, AOB- 
amoA, nirK, nosZ) in regulating №0 emissions are also shown. White squares and 
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white ellipses denote М pools and gaseous М forms, respectively. Black arrows denote 
N transformation directions. The cloud ellipse denotes N2O fluxes. Green ellipses 
denote functional genes that are associated with N transformation processes, while the 
size of the ellipses denote the abundance magnitude of four functional genes found in 
this study (i.e. larger ellipse indicates higher gene abundance). Dark-yellow squares 
denote soil properties that are significantly correlated with gene abundances. Light- 
green lines and orange lines denote significant correlations at the short- and the long- 
term sites, respectively. The ‘-’ sign in the white ellipse indicates negative correlation. 


Insignificant correlations were not shown in this figure. 
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Fig. 4 
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Supplementary information: 


Table 51 Primer information of selected soil functional genes. 


Gene Primers Primer sequence (5'--3') Location Length References 
amoA IF 5 GGGGTTTCTACTGGTGGT 3' 322-249 (Rich et al., 2003; 

AOB-amoA 500bp Levy-Booth et al., 
amoA 2R 5 CCCCTCKGSAAAGCCTTCTITC 3' 802-822 2014) 
CrenamoA 23F 5 ATGGTCTGGCTWAGACG 3' 7-24 : 

AOAIGA 620bp (Levy-Booth et al., 
CrenamoA 616R 5 GCCATCCATCTGTATGTCCA 3' 611-631 2014) 
nirK Е560-589 5'- GGGCATGAACGGCGCGCTCATGGTGCTGCC -3' 560-589 (оуу Boothetal 

nirK 376bp ° 
nirK R906-935 5-ССОСТТССССААСТТСССССТОСТССАСАС -3' 906-935 2014) 
nosZ-F 5'- CGCTGTTCITCGACAGYCAG -3' 1181-1201 

nosZ O0bp (Rich et al., 2003) 
nosZ-R 5- ATGTGCAKIGCRTGGCAGAA -3' 1880-1900 
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Table 52 Reaction programs of quantitative PCR for selected functional genes. 


Gene Primers Thermal cycling conditions References 
1 атоА 1Е : | i 
AOB 95° 30"-40(95° 15"-53° 15"-72° 40") (Rich et al., 2003; Levy 
amoA amoA 2R Booth et al., 2014) 
_ CrenamoA23F Е 
АОА 95° 30"-40(95° 5"-53° 34"-72° 60") (Levy-Booth et al., 
amoA CrenamoA616R 2014) 
nirK F560-589 - 
nirK 95° 30"-40(95° 5"-65° 34"-72° 60") (Levy-Booth et al., 
nirK R906-935 2014) 
nosZ-F 
nosZ 95° 30"-40(95° 5"-56° 34"-72° 40") (Rich et al., 2003) 
nosZ-R 
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Table 53 Soil physiochemical properties (average + standard error in bracket, п = 3) before Ше tracing incubation. SC, SL, SM and SH indicate 


control, low-N, medium-N and high-N treatments in the short-term N addition site, respectively. LC, LL, LM and LH indicate control, low-N, 


medium-N and high-N treatments in the long-term N addition site, respectively. Same lowercase letters following each values indicate 


insignificant differences (p > 0.05) among different N-treatments of the short-term or long-term site. 


Site Treatment N-adding content pH NH4*-N | NOs-N TOC Total N C/N ratio MBC MBN 
kg N hat a! mg kg! gkg! mg kg"! 

SC 0 3.90(0.05)? 4.4(2.2) 7.1(1.7? 25.5(3.3)2 1.9(0.1)* 13.6(0.5)* 300.4(87.8)? 40.0(18.6)* 
SL 35 3.90(0.01)* 6.7(3.5)* 7.5(0.7?  25.8(4.D^ 1.7(0.4)* 15.3(0.6) 138.8(42.8) 20.8(15.8)* 

ааа SM 70 3.84(0.02)^ 6.0(1.2)* 7.0(1.9)*  25.7(1.7? 2.0(0.5)8 13.31.42 136.2(47.3? 17.0(2.9)* 
SH 105 3.80(0.02) 3.0(1.2)* 9.6(0.8? 27.6(2.9)® 1.9(0.3)* 14.7(0.6)* 192.0(14.3)^ 18.5(0.7)* 
LC 0 3.88(0.03)*  5.9(3.5?* 7.1(3.1)*  26.8(2.2)?^ 2.0(0.1)* 13.5(0.3)* 290.7(127.4)* 27.2(17.5) 
LL 50 3.82(0.04) 6.0(3.9)* 8.1(4.5)* 28.8(8.5)^ 2.0(0.6)* 14.3(0.1? 205.6(64.9)® 18.6(4.2)* 

С LM 100 3.71(0.04) 4.4(2.9)* 9.3(1.60 | 35.6(8.2) 2.4(0.4)* 14.7(1.0? 296.0(61.3)* 35.6(12.7)* 
LH 150 3.67(0.07)°  7.5(3.1?* 10.0.8) 31.2(0.2)^ 2.3(0.2)* 13.8(0.7)* 236.8(54.1)8 20.3(3.9)* 
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Table 54 Gross N transformation rates (average + standard deviation in bracket) during the incubation of 8 N-treatments. Kinetics: 0 = zero- 


order, 1= first-order; LSDs or LSD: least significant difference when comparing any two means in the short-term or the long-term site for p = 


0.05. Rates are in mg N kg"! 4". Detailed explanation of the N rates are as follows: Mec: mineralization of recalcitrant organic nitrogen to 


МН.) Гун4-М№ес: immobilization of МНа* to recalcitrant organic nitrogen; Ммиь: mineralization of labile organic nitrogen; Гуна-мав: 


immobilization of NH4* to labile organic nitrogen; Омгес: oxidation of recalcitrant organic nitrogen to МОз Inos: immobilization of NO3 to 


recalcitrant organic №; Охна: oxidation of NH4* to МОз; Dyos: dissimilatory NO3 reduction to NH4*; Аун: adsorption of NH4* on cation 


exchange sites; К мнаа: release of МНД? on cation exchange sites). МЕ Тунчаргоа: net NH4* production; МЕ Туозаргоа: net NO3 production. 


Short-term Long-term 
N rates Kinetics SC SL SM SH LSDs LC LL LM LH LSD. 
Мес 0 1.3(0.3)Е-01  0.7(0.3)Е-01 1.100.4)Е-01  1.1(0.3)E-01 0.031 0.6(0.2)Е-01 0.9(0.5)Е-01 1.1(0.4)E-01  1.4(0.7)E-01 0.046 
Inna-Nrec 1 5.8(3.7)Е-03 2.7(1.5)Е-03 6.4(4.4)Е-03  4.2(2.4)E-03 0.003 4.4(1.5)Е-03  3.3(1.0Е-03  6.2(3.6)E-03  1.1(0.5)E-01 0.026 
Ммаь 1 4.6(2.5)Е-02 8.4(2.9)Е-02  1.1(0.DE-01 7.2(3.7)Е-02 0.032 1.3(0.3)E-01 7.4(4.3)Е-02  1.1(0.5)E-01 . 7.0(6.1)E-02 0.043 
[NH4-Nlab 1 7.5(2.7)Е-03 2.8(1.5)E-03 4.5(2.8)E-03  5.7(3.3)E-03 0.002 3.4(0.9)E-03 2.0(0.6E-03)  5.0(3.7)E-03 1.4(0.4)Е-01 0.020 
Омес 0 1.0(0.9)Е-05 1.5(1.0)Е-05  1.2(1.DE-05  1.0(0.6)Е-05  8.6Е-06 2.2(0.7)Е-05  7.2(5.0)Е-06 = 1.9(1.0)Е-05  1.4(0.4)Е-05 6.4Е-06 
Inos 1 3.8(4.9)E-03 2.5(1.8)Е-02 1.8(1.1)Е-02  1.1(0.9)E-02 0.011 3.0(1.6)Е-02 8.0(7.5)Е-03 5.2(1.4)E-02 5.0(4.0)Е-03 0.011 
Очна 1 1.8(0.1)Е-01 1.9(0.2)Е-01 1.6(0.2)Е-01  1.8(0.2)E-01 0.017 1.7(0.0)E-01 1.6(0.2)Е-01 . 1.8(0.2)E-01 8.0(0.8)Е-02 0.014 
Dyo3 1 3.4(1.6)E-03  1.6(0.DE-02  1.2(0.2)E-02 5.0(1.2)Е-03 0.001 6.0(0.8)E-03 7.0(1.3)Е-03  4.5(2.))E-03  6.9(1.3)E-03 0.001 
Anna 1 1.4(0.5)E-02 5.7(6.1)Е-03 1.1(0.5)Е-02  1.6(0.8)E-02 0.006 1.4(0.4)Е-02  1.2(0.8)E-02  2.6(0.5)E-02 1.0(0.6)Е-02 0.005 
Кунча 1 1.5(0.4)Е-04  5.3(1.9E-05  7.7(7.0)E-05  1.4(1.2)E-04 6.8Е-05  1.1(0.7)E-04 8.4(3.4)Е-05  2.7(0.)E-04 1.3(0.4)Е-04 5.2Е-05 
NET nus prod -2.6(3.0)E-03 -3.1(2.9)E-02 5.4(5.7)Е-02  -1.9(3.4)E-02 0.092 2.2(3.6)Е-02 -3.4(0.4)E-03 8.3(0.6)Е-03 -1.3(0.3)Е-01 0.144 
NET wos. Prod 1.7(0.0)Е-01 1.5(0.0)Е-01  1.3(0.00E-01 1.7(0.1)Е-01 0.029 1.4(0.1)Е-01  1.4(0.DE-01  1.2(0.0Е-01 6.9(0.2)Е-02 0.025 
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Figure legends 


Fig. S1. ЭМ tracing model (Müller et al., 2007) (NH4*: ammonium; NHa*aa; : 
adsorbed NH4*; Nia»: labile organic N; NOs: nitrate; Nrec: recalcitrant organic М. 
Ммес: mineralization of recalcitrant organic nitrogen to NH4*; INH4-Nrec: 
immobilization of МНД? to recalcitrant organic nitrogen; Ммав: mineralization of 
labile organic nitrogen; Гуна-мађ: immobilization of МНД to labile organic nitrogen; 
Омес: oxidation of recalcitrant organic nitrogen to NO; [моз: immobilization of МОз 
to recalcitrant organic М; Oyna: oxidation of NH4* to КОз ; Duos: dissimilatory МОз 
reduction to NH4*; Амна: adsorption of МНА? on cation exchange sites; Кунча: release 


of МН? on cation exchange sites). 


Fig. S2. The measured (scatters, average + standard deviation, n = 3) and modeled 
(lines) N concentration (а, b, d, e) and М enrichment (с, f) in inorganic-N of the 
eight N treatments during the incubation. Symbols (points and lines) in green and 
orange indicate the short-term and long-term N deposition treatments, respectively. 
SC, SL, SM and SH indicate control, low-N, medium-N and high-N treatments in the 
short-term site, respectively. LC, LL, LM and LH indicate control, low-N, medium-N 


and high-N treatments in the long-term site, respectively. 


Fig. S3. Pearson correlation analysis results of soil properties, N20 fluxes, N20 
emission proportions from ammonium (!№О/' МН?) and nitrate pools 


(PN20/PNO»), gene abundances, denitrification/nitrification ((nirK + nosZ)( AOA- 
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атоА + AOB-amoA) ) апа nosZ/nirK gene ratios in Ше short-term (a) and long-term 
sites (b). N204, N2On and №Оа indicate different N2O fractions from autotrophic 
nitrification (a), heterotrophic nitrification (h) and denitrification (d) pathways, 
respectively. Squares in red, blue and yellow (light and medium) indicate NH4*-N, 
NOs-N and organic-N pools (including PN labelled components coming from 
original isotopic addition solutions or from N transformation). Ovals in dark yellow 
and light green indicate soil properties and gene indicators (gene abundance and 
relative ratio), respectively. Overlaps in adjacent indicators indicate significant 
correlations, while insignificant correlations were not shown. There is no meanings 
for the size of the symbols mentioned in this figure. АП the correlations mentioned in 


this figure are positively significant unless otherwise stated. 
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Fig. S3 
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